Nonlinear dynamics of a cigar-shaped Bose-Einstein condensate coupled with a single 

cavity mode 
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We investigate the nonlinear dynamics of a combined system which is composed of a cigar-shaped 
Bose-Einstein condensate and an optical cavity. The two sides couple dispersively. This system is 
characterized by its nonlinearity: after integrating out the freedom of the cavity mode, the potential 
felt by the condensate depends on the condensate itself. We develop a discrete-mode approximation 
for the condensate. Based on this approximation, we map out the steady configurations of the 
system. It is found that due to the nonlinearity of the system, the nonlinear levels of the system 
can fold up in some parameter regimes. That will lead to the breakdown of adiabaticity. Analysis 
of the dynamical stability of the steady states indicates that the same level structure also results in 
optical bistability. 

PACS numbers: 37.10.Jk, 37.10.Vz, 42.50.Pq, 42.65.Pc 



I. INTRODUCTION 

Recently, a lot of investigations have been devoted to 
the combination of ultracold atom physics and cavity 
quantum electrodynamics. Experimentally, the efforts 
culminate in the successful coupling of a Bose-Einstein 
condensate (BEC), that is, a single matter-wave field 
mode, to a single cavity mode in a high finesse opti- 
cal cavity [l], [J Though of a very short history, this 
type of combined system has demonstrated many inter- 
esting phenomena, such as vacuum Rabi splitting 
^/N scaling of the atom-photon coupling [l|, [3], opti- 
cal bistability [1, 01, cavity enhanced superradiance of 
a BEC [5], and most surprisingly, a map between the 
atom ensemble-cavity system and the canonical opto- 
mechanical system is found Besides playing an 
active role in the dynamics, the cavity can also be used 
to characterize the properties of the matter-wave field. It 
has been experimentally implemented to study the statis- 
tics of an atom laser by detecting single atoms falling 
through a cavity [s*]. Based on similar principles, a pro- 
posal to probe the superfluidity-Mott insulator transition 
is brought up [9]. 

One of the most remarkable characteristics of the atom 
ensemble-cavity system is its intrinsic nonlinearity. If we 
are only interested in the mechanical motion of the atoms 
and thus confine ourself to the large detuning limit, the 
atom-photon interaction is of a dispersive nature. The 
atom-photon interaction provides a potential (an opti- 
cal lattice) to the atoms, and meanwhile shifts the cav- 
ity mode frequency. By shifting the cavity mode fre- 
quency, the atom ensemble can influence or even de- 
termine the potential acting on itself. This entails a 
nonlinearity substantially different than the usual atom- 
atom interactions. It is this nonlinearity that lies at 
the heart of previous experimental and theoretical works 
[3. [gI. [tI. [ToI. i3| ■ in particular, the dispersive optical 
bistability (alU"^ 

Along this line, we consider the nonlinear dynamics 



of a cigar-shaped Bose-Einstein condensate in an opti- 
cal cavity in this work. Experimentally, this has been 
achieved and investigated by Esslinger et al. and 
in the weak atom-atom interaction and weak excitation 
limit (or non-depletion limit), it is found that the con- 
densate can be described as a harmonic oscillator and the 
condensate-cavity system maps onto the generic cavity 
opto-mechanical system perfectly. In this work we ap- 
proach the same problem from the perspective of diffrac- 
tion of matter-wave by an optical lattice. We will not 
restrict to the weak excitation limit. We will develop a 
discrete-mode approximation (DMA) for the condensate, 
and by comparing with the full description of Gross- 
Pitaevskii (GP) equation, we find that the finite-mode 
approximation is quantitatively good in the case of weak 
atom-atom interaction. The DMA also enables us to find 
out the steady states of the system, analyze their stabil- 
ity, and elucidate the possibility of breakdown of adia- 
baticity and optical bistability in this system. 

This paper is organized as follows. The system and 
the general formalism are described in Sec. [Tll Then in 
Sec, nil Al we compare the DMA with the full GP dynam- 
ics so as to validate the DMA. In Sec. IIIIBl we address 
the problem of adiabatic evolution of the system, which 
is relevant to optical bistability. Our results are summa- 
rized in Sec. IIVI 



II. GENERAL FORMALISM 

We assume that a cigar-shaped Bose-Einstein conden- 
sate is located inside an optical cavity with its axial 
direction parallel to that of the cavity. The internal 
transition frequency of the two-level modeled atoms is 
LUa, while the relevant cavity mode interacting with the 
atoms is of frequency Uc- The cavity mode is coherently 
driven by a pump laser with frequency ujp and a possi- 
bly time-dependent amplitude r]{t). The atom-pump and 
cavity-pump detunings are denoted as — uja~(-Up and 
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Ac = Wc — Wp, respectively. In the large detuning limit 
and by transferring into the rotating frame at the pump 
frequency, we get the Hamiltonian for the condensate- 
cavity system {h~l throughout in this paper) Q, 

/ I (f 

H= I dx-^Ux)[ — + Vext{x) + UQa^acos^{kx) 

\ 2m d'^x 



+^¥{x)i!{x)ji>{x) 

-t-Aca^a + ri{t){d + a)) + H^. 



(1) 



Here, and are the creation operators for the atoms 
(with mass m) and the cavity photons, respectively. 
Vext{x) — imwyo;^ is the harmonic potential in the axial 
direction with characteristic frequency uj\\ , and giD is the 
effective atom-atom interaction strength in a transversely 
tight-confining trap plj | . Besides the harmonic potential 
Vext (which may also owe to atom-photon interaction, 
as in an optical trap), the atom-cavity photon interac- 
tion provides an additional potential UqO^ acos^{kx) for 
the atoms. The difference is that, the former is static, 
while the latter, being proportional to the intra c avity 
photon number, may be dynamical and quantized [22l |. 
Here [/o — ~9o/^a. is the maximal light shift per pho- 
ton an atom may experience (at an antinode), with go 
being the atom-photon coupling constant. Note that we 
assume that the transverse radius of the condensate is 
much smaller than the cavity mode waist, and therefore 
it is legitimate to neglect the variation of go in the trans- 
verse direction. The cavity mode function in the axial 
direction is cos(fcx), with the wave vector k = 2tt/X. Fi- 
nally, Hf^ accounts for cavity photon decay with a rate 

K. 

In the mean field approximation, we take ^(x,i) ~ 
^(x,i) and d ^ a, i.e., both the matter-wave field and 
the electromagnetic field are described as classical fields. 
The Gross-Pitaevskii (GP) equation for the condensate 



1 d"^ 

- —-j:^ + Vext{x) + Uo\a\^ cos2(A:a;) 



9«'(x,i) 

+giD\'f{x,t)\'j'i>{x,t), 

and the equation of motion for the cavity field is 

da 
'dt 



(2) 



i[Ac + Uo y da;|*(x,t)rcos^(fca;) ) a 
Ka + Tjit). 



(3) 



We can numerically integrate these coupled equations 
and so as to study the dynamics of the combined system. 
To this end, technically, we had better convert the orig- 
inal equations into their dimensionless form. The char- 
acteristic length scale and energy scale of the system are 
f — A/27r and Ur — fc^/2m, i.e., the period of the optical 
lattice divided by tt and the atomic recoil energy, respec- 
tively. We thus rescale position and time as x = and 



t = t/ujr respectively, and accordingly the condensate 
wavefunction ^'(a;,i) = N/^^{x,i), so that the scaled 
wavefunction is normalized to unity, / di|\E'(i,t)p — 1. 
Here N is the total atom number. Equations ^ and ([3]) 
then convert to 

= (-i.+.^i^ + Uo\a\^cos^S: 
dt \ d^x 



+5i*(i,t)n*(i,i), 



(4) 



-Z = -^i(Ac + Nijo dx\^{x,i)\'^ cos^ x)a 
dt J 

— ka + ijlt). 



(5) 



Here we introduced the rescaled dimensionless quantities 
(Uo, Ac, K.,ri) = {Uq, Ac, K,r])/ujr, the rescaled harmonic 
frequency lu = y^mcij^^^/2ajr, and the atom-atom inter- 
action strength g ~ Ngiu/oJrS,- 

At this point, a further simplification is possible. Note 
that experimentally the cavity damping is much faster 
than the mechanical motion of the condensate. The for- 
mer occurs at a rate of n, while the latter, roughly speak- 
ing, is on the order of LUr (see the discrete-mode approx- 
imation below) . In the experiment of Esslinger et al. 0] , 
K = 27rxl.3MHz, and Ur is around 27rx3.75 kHz (^"^Rb), 
so the cavity decay is almost three orders of magnitude 
faster than than the condensate motion. We thus can 
safely assume that the cavity field follows the condensate 
adiabatically and solve the photon number as 



riph = |a(i)P = 



+ {Ac + NUoicos^ i))2' 



(6) 



where {cos^ x) is defined as / di|\P(i, i)p cos^ i, and has 
the meaning of the overlap between the cavity mode in- 
tensity and the condensate density distribution. Obvi- 
ously, it is bounded from above and below, < (cos^ x) < 
1. Equation ^ gives the dependence of the intra cav- 
ity photon number, i.e., the intra cavity optical lat- 
tice depth on the condensate. The cavity-pump detun- 
ing is shifted by the condensate to an effective value 
Acff = Ac + NUo (cos^ x) . 

Substituting Eq. ^ into Eq. (g]) , we get a GP equation 
for the condensate, the most prominent feature of which 
is that the potential acting on the condensate depends in 
a highly non-local and non-linear way on the condensate 
itself. This is to be compared with the non-linear term in 
the usual GP equation, the atom-atom interaction, which 
is nonlinear but local. If this cavity induced nonlinearity 
is to exhibit itself apparently, the influence of the conden- 
sate on the cavity field should be strong. To this end, it 
is desirable to have k NUo (with an appropriate Ac). 
Only under this condition, the condensate is able to shift 
the cavity into or out of resonance with the pump [see 
Eq. ©J and so as to influence the cavity field drastically. 
It is worth noting that this condition has been fulfilled 
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with both Fabry-Perot cavities [1, H, 13 ^^'^ '"i^S cavities 

Although a brutal integration of Eq. ([4|) is possible, 
and that captures all the possible atomic modes simul- 
taneously, we would like to introduce the discrete-mode 
approximation (DMA). The advantage of DMA is that 
it captures the physics and is technically simpler than 
the full GP equation. To explain the basic idea, let us 
consider the extremal case to = and g = 0, i.e., a non- 
interacting homogeneous condensate. The ground state 
of the condensate is the zero momentum state \p — 0). 
The effect of the intra cavity optical lattice is to diffract 
this state into the state (|p = 2k) + \p ~ — 2fc))/V2 
(diffraction into the asymmetrical superposition state 
is forbidden by the parity symmetry of the potential), 
which can then be further diffracted into {\p = 4fc) -|- |j5 = 
— 4fc))/-\/2 (or back into \p = 0)), and so on. Thus the 
dynamics of the system actually confines to some discrete 
atomic modes. This stimulates us to introduce the DMA. 

In the actual case, uj and 5^0, the ground state 
\['g(a;) is not strictly the zero momentum state. However, 
usually the length of the condensate is much larger than 
the period of the optical lattice, thus by the long range 
coherence of the condensate and the uncertainty relation, 
we infer that the momentum distribution of the conden- 
sate is well localized around the origin with a width much 
smaller than 2k. Thus the picture for the non-interacting 
homogeneous BEG carries over to the actual long 
as in the period we are concerned with, the atomic diffrac- 
tion process overwhelms other effects caused by the finite 
trap potential and atom-atom interaction. We can take 
into account the small but finite momentum spread of the 
initial state by defining the discrete modes as — 'i'g {x) , 
= V2cos(2A:a;)*g(x), 02 = V2cos{4:kx)-^g{x), and 
generally for n G N, (pn — V2 cos{2nkx)'^ g{x) In 
the momentum representation, (/>„ is a coherent superpo- 
sition of two peaks localized at ±2nfc respectively. For 
the n-th mode, we introduce the creation (annihilation) 
operator cjj (c„), with common boson commutation rela- 
tions. Note that these mode functions are not strictly or- 
thonormal, so this assignment of commutation relations 
is a bit problematic. However, as long as the ground state 
^'g(a:) varies slowly on a length scale of A/2, the error is 
negligible. 

The DMA assumes that the condensate dynamics is 
confined in the subspace spanned by {4>n}- Numerically, 
we find that it is enough to cut off at 02 • Higher modes 
play a minor role due to the quadratically increasing 
eigenenergies. Substituting '^{x) = X]^=o '^*'?^*(^) 
Eq. ([T]), and neglecting the atom-atom interaction, we 
get the Hamiltonian in the DMA 



H — 4:UJrc[ci + IGuJrcl^'^ + ^fl^fl (\/2(cJci + c[cf)) 

+c\c2 + 4ci -I- 2A^) + Acfltfl + r]{a + a^) + i?(7) 



In our calculations, we have used the equations: 

dx4>l ( "2^^ + ^ext 1 00 = eo, (8a) 

j dx(t>*„^ 2m ^ ^ext^ 0n ^ 5mn{St) + An^UJr), (8b) 

dx(t>*^ COS^ {kx)4)n ^ ^{5m.n-l + ^m,n+l), (8c) 

and dropped the constant term proportional to ep- These 
equations become exact in the non-interacting homoge- 
nous case. In the DMA Hamiltonian, the atom-photon 
interaction can be interpreted in two ways. For the 
atoms, they are diffracted between adjacent modes with 
coupling strengths proportional to the photon number; 
For the photons, their frequency is renormalized propor- 
tional to the interference of the atomic modes. 

The atom-atom interaction we neglect contains terms 
like cjcicjco which shift the energies of the ex- 

cited modes, and more importantly, it mixes the dis- 
crete modes we retain and other modes we discard in a 
four-wave mixing way. Therefore, the atom-atom inter- 
action impairs the discrete-mode approximation. How- 
ever, since the energy scale of the atom-atom interaction 
is on the order of the chemical potential of the conden- 
sate /X = (3(7^/4)^/^ (in units of uj^, in the Thomas- Fermi 
regime) , it may be legitimate to neglect it in a time scale 
of if M 4, the free frequency of the first excited 
mode in units of lOj.. In Esslinger et al's experiment, 
Aojr ~ 27r X 15 kHz, while the chemical potential of the 
condensate is about 27r x 2.4 kHz. 

From the Hamiltonian ([7]), by taking the mean field 
approximation Ci ^ \[NZi^ a ^ a, and integrating out 
the cavity field, we get the equation of motion for the 
condensate 



t Z = H{nph)Z = [Hi 
at 



'l'phH2) Z, 



(9) 



with Z = {Zq, Zi, Z^)^ ^ and 



Hi 



16 




(10) 



The photon number is 



k2 + (A^ + ANZ^H^Zf ' 



(11) 



where = Ac + NUq/2 is the rescaled effective cavity- 
pump detuning with the condensate in its ground state 
[Z = (1, 0, 0)^]. We would also fike to define the general 
effective cavity-pump detuning Ae// = A^ -I- AN Z"^ H2Z . 
Equation (fTT|) is just Eq. ^ in the discrete-mode ap- 
proximation. 
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III. APPLICATION OF THE FORMALISM 

In this section, we apply the formalism developed in 
the preceding section to two different cases. The first 
case deals with a constant pump, which acts on the BEC- 
cavity system abruptly, the evolution of the BEC-cavity 
system is studied. The second case considers the scenario 
that the pump is turned on slowly and smoothly. We aim 
to study the evolution of self-sustained steady state of the 
BEC-cavity system as the parameters varied. 



A. Constant pump: validity of the discrete-mode 
approximation 

Assume that initially the BEC is in its ground state, 
and the cavity is the vacuum. Then at i = a beam 
of laser with constant frequency and constant amplitude 
is projected onto the cavity. In a time scale of the 
intra cavity field and hence the intra cavity optical lat- 
tice builds up, which then diffracts the BEC into higher 
momentum states. In the position representation, the 
condensate is deformed by the optical lattice potential. 
The point is that, as the condensate deforms, its over- 
lap with the cavity mode and hence the effective cavity- 
pump detuning varies, which in turn leads to the vari- 
ation of the lattice potential. Though atom diffraction 
by an optical lattice has been extensively studied both 
theoretically and experimentally since two decades ago 
0, m, [13, [Tsl . [Toj , in all previous works the optical lat- 
tice is either static or just a time-dependent parameter. 
In contrast, here the lattice potential has to be treated 
as a dynamical variable. 

We have simulated the dynamics of the BEC-cavity 
system with both the CP equation and the discrete- 
mode approximation. The CP simulations are shown in 
Figs. [Ija)-(c), with oj ~ 0.01 being fixed, and g — 100, 
50, and 10, respectively. In each figure, we present the 
evolution of the photon number Uph and the normalized 
populations P„ = p on different atomic modes. 

In one dimension, the Thomas-Fermi approximation is 
valid if 5 ^ \J2txuj With this criterion, we find that 
the three cases are all deep in the Thomas-Fermi regime. 
Thus, the atom-atom interaction plays an important role 
in determining the ground state of the condensate. In 
particular, the radius of the condensate is enlarged from 
(l/cj)i/2 to Rtf = (4g/3w)i/3 (-^^ ^^its of ^ = A/27r). 
For all the three cases, we have Rtf ^ ti", therefore the 
basic condition for DMA is satisfied. 

The simulation based on DMA is presented in 
Fig. [Ijd). Note that in DMA, the atom-atom interac- 
tion and the initial ground state are irrelevant. Thus, 
Fig. [IJd) is a common approximation to Figs, ma)-(c). 
By comparing Fig. [Ijd) with Figs. [IJa)-(c), we see that 
to a time scale as long as ijJrt — 20, the discrete-mode ap- 
proximation agrees with the full GP equation very well. 
Moreover, the agreement improves as the atom-atom in- 
teraction decreases. The DMA recovers the 5 = 10 case 



[Fig. [Ijc)] almost exactly. By examining the difference 
of Figs. [Ija)-(c) with Fig. [Ijd), we identify two conse- 
quences of the atom-atom interaction. Firstly, we see in 
Figs.[IJa) and[IJb) that as time goes on, the total pop- 
ulations on the discrete modes we take decreases, and 
the time scale of this process increases with decreasing 
atom-atom interaction. Therefore, it is reasonable to at- 
tribute this phenomenon to the atom-atom interaction, 
which populates the atoms outside the modes we take. 
Secondly, the atom-atom interaction shifts up the fre- 
quencies of the excited atomic modes, and thus shortens 
the oscillation period of the system 21 1 . 

To illustrate the second point, we plot in Fig. O the 
evolution of the effective detuning Ag// in Figs.[Ua)-(d). 
The oscillation of Ae/ / is directly related to that of the 
photon number. We see clearly that the periods in GP 
simulations are a bit shorter than that in DMA, and the 
stronger the atom-atom interaction, the more apparent 
this effect. Figure [5] also helps us understand the double- 
peak structure of the photon number line in Fig. [1] Each 
time Ae// crosses zero, rtp/i attains its maximal value 
f]^ Ik} . 



B. Slowly varying pump: breakdown of 
adiabaticity and optical bistability 

In this subsection, we proceed to consider the scenario 
that initially the BEC is in its ground state and the cavity 
is in the vacuum, and subsequently the pump is slowly 
and smoothly turned on. In contrast to the situation 
above where the system is exposed suddenly to a finite 
pump and experience oscillations indefinitely, here we 
are concerned with the adiabatical following of the self- 
sustained steady states. By a self-sustained steady state, 
we mean that the two sides are stationary and consistent 
with each other: the condensate is in the ground state de- 
termined by the optical lattice (plus the harmonic trap), 
and meanwhile the optical lattice (its depth) is deter- 
mined by the condensate configuration, i.e., fulfilling ([B]) 
or pT|) . Naively, it is expected that as long as the pump is 
ramped up sufficiently slow, the BEC-cavity system will 
follow the steady states all the way. However, now we are 
dealing with a nonlinear system and the nonlinearity may 
give rise to the breakdown of adiabaticity [23l.[2^.[25l.[2g|. 
That is what we are interested in. 

For our purpose, the key issue is to map out all the pos- 
sible steady states as the pump strength varies. Although 
this is a bit difficult (time-consuming numerically) in the 
GP formalism, it can be done very conveniently with 
DMA. In the framework of DMA, a steady state cor- 
responds to a solution of Eq. ([9]) such that 

Zo(i) = Zoe-'^«*. (12) 

Or equivalently, 

H{nph)Zii ^ [Hi + nphUi) ^ EqZi^, (13) 
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FIG. 1: (Color online) Simulations of the BEC-cavity dynamics, (a), (b) and (c): with Gross-Pitaevskii equation; (d): with 
discrete-mode approximation. In each panel, the black solid line is for the photon number, while blue dashed line, red dotted 
line for the normalized populations on the atomic modes (j>o and (pi , respectively (populations on the third mode are too small 
and not shown). The parameters are iV = 4.8 x lO", Uo = 0.25, {R,A'^,ij) = (0.4,0.4,2.4) x 10^ uj = 0.01, and g = 100, 50, 
and 10 in (a), (b) and (c), respectively. 
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FIG. 2: (Color online) Evolution of the effective detuning 
Ae//. The dashed line, dashed dotted line, dotted line, and 
the solid line correspond to Fig. [TJa)-(d), respectively. 



which defines a nonlinear eigenvalue problem because the 
Haniihonian H depends on the eigenstate Zq through 
Hph- It is impossible to invoke the concepts and tools 
in linear algebra to solve this problem. Our strategy is 
to first take an arbitrary trial photon number ritj., solve 
the ground state of the Hamiltonian H{ntr), Ztr, and 
then substitute it into Eq. (fTT|) to get an output pho- 
ton number Hout- If nout = ntr, then the solution is 
self-consistent and a steady state is obtained. Note that 
Uout, as a function ointr, is continuous and moreover, it is 
bounded both form above and below, < Hout < fj^/K^- 
Therefore, if we scan ritr from to 77^ /k^, by the princi- 
ple of continuity, there must be at least one point where 
nout = n-tr- This guarantees the existence of steady state 
solutions (clearly the same strategy and arguments apply 



also in the GP formalism). Furthermore, because of the 
nonlinearity of the system, the possibility of more than 
one steady state solutions for a given set of parameters is 
expected. As observed and argued in Refs. [23,, 2^, 
the existence of more eigenvectors than the dimensional 
of the Hilbert space is unique for nonlinear systems. 

In Fig.[3]Ja), we depict the output photon number Uout 
as a function of the trial photon number ntr as the pump 
strength fj varies while other parameters are held at con- 
stant. For each curve, each of its intersections with the 
dotted line Uout — ntr corresponds to a steady state. 
We see that for fj < 960 and fj > 1440, there is only 
one steady state solution, while in an intermediate range, 
there can be three (for example, fj = 1200). This fact is 
made more evident in Figs.[3]^b) and[31[c), where we plot 
the photon number at steady state Ust and the nonlinear 
eigenvalue Eq as a function of fj directly. In the inter- 
val {fj2,fji) = (1016,1296), each fj corresponds to three 
steady states and hence three rigt's and three E^s. The 
folded level structure in Fig. [3]Jc) is similar to the looped 
level structure in Refs. [23, 25j- 

With these figures as maps, we return to the scenario 
again. As the pump is slowly ramped up from zero, the 
BEC-cavity system follows the curve OH I in Fig. [3fb) 
adiabatically. Things go on like this until the pump ex- 
ceeds 771, where the curve terminates. Interpreted in 
terms of Fig. [2Ia), the system first evolves from A to 
B to C, then at some point C coincides with D and then 
they disappear all together. Similar interpretation can be 
made in terms of Fig. Efc). As 77 crosses 771, the system 
cannot jump discontinuously from / to K . Therefore, be- 
yond this critical value of pump, there is no way for the 
BEC-cavity system to evolve adiabatically anymore, no 
matter how slow the external pump is varied. We thus 
expect large amplitude oscillations to appear abruptly as 
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FIG. 3; (Color online) (a) The output photon number Uout as a function of the trial photon number ntr- The intersection 
points {A — G) of the curves with the diagonal dotted line correspond to steady states. The parameters are A'' = 4.8 x 10*, 
Uo = 0.25, (k, A^) = (0.4, 1.2) x 10^, and from bottom to up ?? = (0.72,0.96, 1.20, 1.44, 1.68) x 10^ respectively, (b) Photon 
number at steady state rist and (c) nonlinear eigenvalue Eo as a function of the pump strength rj. Parameters other than fj are 
the same as in (a). The coordinates of I and J are (»7i,nsti) — (1296,2.416) and (f}2,nst2) ~ (1016,6.36), respectively. 



in Fig. [T] This is verified with simulations both based on 
the GP equation and DMA, as shown in Fig. [H There 
we assume that the pump is ramped up and down with 
a Gaussian profile. Both the simulations reveal clearly a 
two-stage feature of the history of the BEC-cavity sys- 
tem. Before a critical moment tf. ~ 40 (corresponding to 
Vc ~ fji), the photon number varies smoothly and coin- 
cides with the prediction based on DMA perfectly. This 
indicates the system evolves adiabatically along the line 
OHI. However, just across this point, pronounced os- 
cillations with large amplitudes and small periods occur. 
Note that here the breakdown of adiabaticity is related 
to the disappearance of a nonlinear eigenstate. This is 
to be compared with its counterpart in linear systems, 
which usually results from level-crossing. While the lat- 
ter can be avoided in principle, the former cannot. In the 
viewpoint of state preparation, to prepare the system on 
the line KL, we cannot take the route we take here in 
the parameter space. 

The structure of the curve in Fig. reminds us of 
optical bistability [2^. This motivates us to study the 
dynamical stability of the steady states. We would like 
to stress that adiabaticity and stability are two different 
but closely related issues. Only those steady states which 
are dynamically stable can be adiabatically evolved [1^ . 
Let Z{t) — Zo{t) + 6Z{t), with 6Z being an infinitesimal 
derivation from the solution Zq. Substituting this into 
the differential equation ([9]), and making use of the fact 
i-^Zo = [Hi + H2nph{Zo, Zq)) Zq, we get to the first 
order in 5Z, 



i—dZi = 
dt 
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FIG. 4: (Color online) Evolution of the cavity photon num- 
ber Uph as the pump is varied as fj{t) = fimax exp{—{ujrt — 
50)^/(7^), with the peak amplitude fimax = 1440 and the pulse 
width a = 31.25, and other parameters being the same as 
in Fig. [3l (a) simulation with GP equation {uj = 0.01 and 
<; = 10 as in Fig. [Tic)); (b) simulation with DMA, the dotted 
line dictates the photon number of steady states drawn from 
Fig. [S^b). Note the sharp edges around ujrt = 40. 



Here the subscript means taking value at [Zq, Zq). As- 
suming 

SZ{t) = e-'^o* (ue-'"* - v*e'^') , (15) 
and substituting this into the equation above, we get 



(14) 



(16) 



7 



where 

B -A-B + Eo )■ 

Here the 3x3 matrices A and B are defined as 



A = Hi+nph{Zo,ZQ)H2, 




0<i,j,k<2. (18) 



The dynamical stability of the steady states is deter- 
mined by the eigenvalues of the matrix M. If all the 
eigenvalues are real, the steady state is dynamically sta- 
ble; otherwise, it is dynamically unstable. 

In Fig. El we examine the stability of the steady states 
on the curve OHIJKL by plotting the maximal imag- 
inary part Im = maxi{/ra(Ai)} of the eigenvalues of 
M (as a real matrix, its eigenvalues appear in conju- 
gate pairs). The curve is parameterized with the photon 
number Ust- We see that for Ust < risti {nst > nst2), 
which corresponds to OHI (JKL), the eigenvalues are 
all real and thus steady states on these segments are sta- 
ble. While for Usti < rist < nst2, which corresponds to 
/J, some eigenvalues are complex and thus steady states 
on this segment are unstable. This is consistent with the 
usual observation that the negative sloped parts are al- 
ways unstable though the positive sloped parts are not 
necessarily stable [23, HI]- As noted previously, stabil- 
ity is a prerequisite of adiabaticity. Here, it is implied 
that we can propagate adiabatically the BEC-cavity sys- 
tem along the line OHI (as we do above), but it is not 
feasible to do that along IJ. 




FIG. 5: Maximal imaginary part Im ~ maxi{/m(Ai)} of 
the eigenvalues of M. The points on the curve OHIJKL 
in Fig. |3][b), which is parameterized with Uat, are examined. 
The interval where Im is nonzero is {nsti,nst2)- 



The analysis of stability confirms the possibility of op- 
tical bistability in the BEC-cavity system. In previous 
experiments dealing with a cloud of thermal cold atom 
gas [5], optical bistability has been discussed and ob- 
served, but an analysis of stability is lacking. If there are 
some dissipation channels, it may be possible to switch 
between the upper and lower stable branches and observe 
the typical hysteresis loop of optical bistability in the 
BEC-cavity system. Here we can only simulate the co- 
herent behavior of the system and thus this phenomenon 
is beyond our means. However, the mechanism underly- 
ing optical bistability still manifests itself in the form of 
breakdown of adiabaticity. 



IV. CONCLUSIONS 

We have investigated the nonlinear dynamics of a 
BEC-cavity system. This system is characterized by 
the nonlinearity that the optical lattice acting on the 
condensate depends on the condensate itself in a highly 
non-linear and non-local way. By borrowing ideas from 
atom optics, we develop the discrete-mode approxima- 
tion (DMA). The philosophy of DMA is that the dynam- 
ics of the condensate confines to some discrete standing- 
wave modes. It is numerically verified that in the limit 
of weak atom-atom interaction, DMA agrees with the 
full CP equation dynamics perfectly. Our emphasis is on 
the adiabatical evolution of the BEC-cavity system. The 
DMA facilitates us to solve the nonlinear eigenstates of 
the system and analyze their stability. Due to the non- 
linearity of the system, the nonlinear energy level can be 
folded in some parameter regimes. This structure gives 
rise to both the breakdown of adiabaticity and optical 
bistability in this system. As an open problem, we con- 
jecture that these two effects are complementary to each 
other. If the dynamics of the system is purely coherent, 
the breakdown of adiabaticity will show up; otherwise, if 
dissipations such as atom loss are strong, optical bista- 
bility may take the place. 
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